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J. Leon and A.V. Mikhailov + 
Physique Mathematique et Theorique, CNRS - UM2, 34095 MONTPELLIER (France) 
+ Landau Institute for Theoretical Physics, MOSCOW (Russia) and 
Dept. Applied Math., Leeds University, LS2 9JT, LEEDS (U.K.) 

Stimulated Raman scattering of a laser pump pulse seeded by a Stokes pulse generically leaves 
a two-level medium initially at rest in an excited state constituted of static solitons and radiation. 
The soliton birth manifests as sudden very large variations of the phase of the output pump pulse. 
This is proved by building the 1ST solution of SRS on the semi-line, which shows moreover that 
initial Stokes phase flips induce Raman spikes in the pump output also for short pulse experiments. 



Introduction. Stimulated Raman scattering (SRS) of 
high energy Laser pulse in two-level medium has been in- 
tensively studied these last years (see e.g. the review 0), 
more especially after the experiments by Driihl, Wen- 
zel and Carlsten Q showing that spikes of pump radi- 
ation (Raman spikes) occur spontaneously in the pump 
depletion zone, see also ||. Consequently nonlinear Ra- 
man amplification was used as a means to observe at 
the macroscopic level the fluctuation of the phase of the 
initial Stokes vacuum 0. 

The SRS equations possess a Lax pair j|] and, as a con- 
sequence, theoretical and experimental works were partly 
devoted to the search of the Raman soliton predicted in 
|^| . In particular it has been believed || that the spike 
of pump radiation observed in Q is a soliton which, in 
the inverse spectral transform scheme (1ST), would be 
related to the discrete part of the spectrum. However it 
has been proved that the observed Raman spike is not 
a soliton but merely a manifestation of the continuous 
spectrum 0. 

Then the question of the creation a genuine Raman 
soliton is still an open problem (for 20 years now). We 
show here that the SRS equations, solved as a boundary 
value problem on the semi-line, do induce the generation 
of solitons by pair, and that, after the passage of the 
pulses, the solitons are static in the medium. Although 
the solitons are not seen in the output pump intensity 
(where only Raman spikes are seen), the signature of 
the soliton birth consists of very large fluctuations of the 
phase of the pump output, which remains to be experi- 
mcntaly observed. 

Last we prove that the Raman spikes, observed in long 
pulse experiments (pulses with fwhm of 75 ns in Q and 
15 ns in ||), where damping is essential in the model, 
must occur also for short pulses (no damping), and give 
the condition for their creation. 

The efficiency of the method of solution serves not only 
as a basis for seeking the Raman soliton, but also to 
evaluate exactly the result of any Raman amplification. 
Indeed the SRS equations are solved for arbitrary val- 
ues of the input (x = 0) pump and Stokes pulses on a 
medium with any initial (t = 0) state. The original prob- 
lem is reduced to a Riccati equation (O) and linear in- 



tegral equations (|12|) . In many interesting cases the Ric- 
cati equation ( [Tl| ) can be explicitly solved (see (|l5|)) and 
important physical information can be obtained without 
solving equations (|l2|), such as the output jl^). 

The SRS equations. For a medium of length L, a de- 
tuning Au and a phase missmatch Ak, in the dimen- 
sionless variables x = Z/L, t = (cT — rjZ)/L, where 
Z and T are the physical variables, the SRS equations 
for short pulses (infinite dephasing time) in the slowly 
varying envelope approximation read M 
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The phase missmatch is represented there by the dimen- 
sionless parameter k = ^L(Ak — 2 Aw). The envelopes 
aT, (frequency ojl) and as (frequency u>s) are scaled to 
the peak intensity I m of the input laser field and the 
dimensionless coupling constant is given by 
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where a' Q is the differential polarizability at equilibrium, 
c/ry is the light velocity in the medium of transition fre- 
quency u>vi and N is the density of oscillators of mass m. 
We have considered in ([!]) the cooperative interaction of 
all fc-components with the medium and then the input 
(x = 0) values of the light waves are also function of k 
sharply distributed around k = 

a L (k,0,t) = J L (k,t), a s (k,0,t) = J s {k,t). (3) 

Finally the electromagnetic field is obtained as 
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The system ([!]) can be solved on the infinte line for 
arbitrary asymptotic boundary values Q and it has been 
proved in M that this allows the interpretation of the 
experiments of H- We will demonstrate that it can be 
solved completely on the semi-line x £ [0, oo), which ac- 
tually furnishes the solution at any point x = L and 
hence solves the finite interval case as a free end prob- 
lem. 
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Lax pair. The time evolution in (|l]) results as the 
compatibility condition Ut — V x + [U, V] — for the Lax 
pair (x > 0, t > 0) 



f x = Uip + ikipa 3 , (4) 
(fit = Vip + w - lk ^ x ^le lk ^ x (5) 

where the dispersion relation fi is x-independent and 
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V^(fc) is discontinuous, therefore only the limits = 
V(k ± iO) arc defined on the real fc-axis. Such is also the 
case for f2, and will have to be determined from the 
boundary conditions on the solution ip. 

Scattering problem. Following standard methods in 
spectral theory, see e.g. |l(|, we dehne the solutions ip + 
and f~ according to 
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These solutions obey the reduction <ft{k) = icr 2 f 2 (k), 
the Riemann-Hilbcrt relations 
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which define the reflection coefficient p and the trans- 
mission coefficient r (p*(k) stands for p(k)). Note that 
det{<p(x)} — det{<p(cx))} = 1, hence for real k |r| 2 = 

The column vectors ip^ and <p 2 are entire functions in 
the fc-plane, with good behaviors as k — > 00 in the upper 
half- plane for ipf (and the lower one for <f 2 ). The vector 
ip 2 is meromorphic in the upper half-plane with a finite 
number N of simple poles k n (and tp^ in the lower one 
with poles k n ). Consequently r and p have meromorphic 
extensions in the upper half-plane where they possess the 
N simple poles k n (the bound states locations). 



Solution of SRS. From the boundary data (||) we have 
readily (set Aw = for simplicity because it can be scaled 
off in as) 
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After careful computations of the boundary values of (J^), 
using separately f + and ip~ , we obtain both limits + 
and f2~ of £1 and the following time evolution of the scat- 
tering coefficient p(k,t): 
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The reconstruction proceeds through the solution of 
the Riemann-Hilbert problem (H) which reads 
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where C+ (resp. C_) is a contour from —00 + iO to 
+00 + iO (resp. from —00 — iO to +00 — iO) passing over 
all poles k n of p{k) (resp. under k n ). Remember the 
notation p*(X) = p(X). 

Since m and 4> are continuous bounded functions of 
k, the coefficients of the Riccati equation ( |TT| ) are ana- 
lytic functions in the upper half fc-plane, where precisely 
p is meromorphic at time t = 0. Consequently p remains 
meromorphic at later times, which allows to prove that 
the reconstructed fields do obey the system of equations 
(|]), or in short that the reconstructed cp does obey (|]). 
Indeed from (|l2|); the Cauchy theorem leads to 
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The formulae (gjj) and fll2j) with ((Hf furnish the solu- 
tion of (Q), in particular the output field values at x — L 
(whatever be q for x > L). In the limit L — > 00, and by 
(^), the solution becomes explicit, and for instance the 
output pump reads 
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Evolution of laser pulses. We are interested here in 
the physically relevant case when the medium is initially 
in the groud state, that is p(k, 0) = 0. Next we choose 
a Stokes wave as a portion of the pump wave, both with 
Loretzian lineshape, namely, for k > 0, 
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n(k 2 + n 2 ) ' 



J s = e 



where A{t) is the pulse shape of duration t m (t > t m => 
A(t) = 0). The parameter 7 (real positive constant) mea- 
sures the ratio pump/Stokes inputs. Considering the case 
when the input Stokes wave experiences one phase flip at 
some time to, we choose the phase 6{t) — for t < to and 
0(t) = vr for t > t . 

It follows from contour integration that the Riccati 
equation (uW can be rewritten as 
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Taking into account the phase flip, the general solution 
with zero initial datum reads 
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cosh((5 — So — 7) — po sinh(<5 — So) ' 
with po — p(k,to), So — S(k,to) and the definitions 
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Soliton generation and signature. Let us consider first 
the case with no phase flip, that is to > t m . The spectral 
transform p, first function in (|l5|), has an infinite set of 
moving poles k n (t), n 6 Z, given by 
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which are associated with solitons as soon as they lie in 
the upper half plane (in the lower half plane they are the 
resonances). 

As t evolves, and for a given linewidth n, these poles 
move from the point —in, and they may cross the real 
axis (and generate solitons) if Tit) is large enough, which 
means enough energy in the input pulses. Moreover, since 
k n = —k- n —i, solitons are created by pair. After the 
passage of the pulse, (t > t m ), m and <j> vanish in (11), 
T(t) is constant, and hence the whole solution becomes 
t- independent. Consequently the laser pulses leave in the 
medium a finite number of static bi-solitons. 

The question of their observation is not straightfor- 
ward. Indeed, taking for simplicity the case of an infi- 
nite medium, the pump output is given in ( |l4| ) where 
the coefficients p and r are understood for real values 



of k (the essential missmatch). The main point here is 
that both 1/t and pjr are holomorphic functions ion the 
upper half-plane of fc, and hence an integrated intensity 
like / dk\aj\ 2 will not show the presence of the poles k n . 
Moreover, if at t — t s a resonance crosses the real axis at 
k s to become a bound state (soliton birth), the output 
intensity at that particular value k s will simply show a 
full depletion (l/r(fc s ) = and |p(fe s )/r(fc s )| = 1). 

However, the phase of the transmission coefficient r(fc) 
diverges for all real k precisely at t — t s . Indeed the trace 
formula (with no pole on the real axis i.e. just before pole 
crossing) reads for k 6 R prf 
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where the integral over A is understood as the Cauchy 
principal value, and where the product concerns only the 
bound states if any. It is then clear that, for every k k s , 
the phase of r(k) diverges as t — * t s because p(A) has a 
pole in A = k s . The observation of true Raman solitons 
should go then with a measure of the phase of the pump 
output. 

We remark here that, in the sharp line limit k — > 
and with absense of a flip in the boundary condition, the 
infinite set of poles k n lie in the upper half-plane and 
\k n \ is proportional to T(t) (|l7j), which implies that the 
solution becomes self-similar (the self-similar solution of 
SRS is analized in fjj|). 

Stokes phase flip and Raman spike. As shown in , 
the Raman spike observed in Q and largely studied later, 
occurs at time t r for which p(k,t r ) = 0. Indeed, then 
|r(&;,£ r )| = 1 and from (|lj) the pump (intensity) is fully 
repleted. 

The second expression in ( |l5| ) gives that p vanishes at 
t r for 
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which relates to (instant of phase flip) and t r (instant of 
Raman spike) to the pulse area through (|l|). 

To get a simple insight in that formula, we may con- 
sider the dominant term k — (resonance) for which <5 
and S r — So are real real positive. Then a Raman spike 
will be allowed if (at k = 0) 

coth(<5 m — S ) > coth(£ ) — 2tanh7 > 0, 

where 5 m — S(k,t m ), and t m is the pulse duration. This 
is obtained without damping, which suggests the search 
of Raman spikes in short pulse experiments (short with 
respect to the medium dephasing time). 

We note finally that the occurence of a phase flip alters 
the denominator of p, which can be used to change the 
discrete spectrum and hence to modify the soliton part 
of the solution. 
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Conclusion. The following results have been ob- 
tained: 

1 - The proof that genuine Raman solitons do occur in 
SRS experiments and cannot be detected in the energy 
of the output but rather in its phase divergence. 

2 - The spontaneous generation of solitons out of vacuum 
in an integrable system. 

3 - The complete solution by 1ST of the boundary value 
problem for SRS on the semi-line. This solution relies on 
the analytical properties of the nonlinear Fourier trans- 
form p(k, t) which are conserved in the time evolution. 

4 - The proof that Raman spikes are generic as they oc- 
cur as well for short duration pulses, as suspected in the 
experimental work 

5 - The proof that in the sharp line limit and with the 
proportional pump/Stokes input on a medium initially 
at rest, the solution is self-similar. 

Comments. 1 - The 1ST solution of SRS on the semi- 
line has been considered in jl3| for the sharp line case. 
Following the same procedure in our inhomogeneously 
broadened case, and in particular neglecting some terms 
in the spectral data evolution, would lead to an evolution 
for p with non-analytic coefficients, in contradiction with 
the necessary meromorphy of p. 

2 - 1ST on the semi-line for the nonlinear Schrodinger 
equation has been considered in [jlJJ , and the problem is 
considerably more difficult mainly because of the neces- 
sity of the knowledge of additional constraints in x = 
(e.g. the time derivative of the field). 

3 - Boundary conditions for nonlinear evolution equations 
are known to be a source for soliton generation Jl5| [63 , 
but this is the first instance of an integrable system which 
shows explicitely the spontaneous generation of solitons, 
as motion of poles from the lower to the upper half-plane, 
from resonances to bound states. 

4 - At the moment we are not able to extract, from the 
explicit form of the spectral data, the explicit form of 
the potential (the medium excitation field) because on the 
semi-line no pure soliton exist (would p vanish on the real 
axis that it would vanish everywhere and no bound states 
would be allowed). However the method provides the ex- 
plicit form of the physically relevant quantities which are 
the light pulses output values. 

5 - The manifestation of the soliton birth as sudden large 
variation of the phase of the output was quite unsus- 
pected and should lead to reconsider the measures of the 
output in SRS experiments. 
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